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The insulating ground state of 5d transition metal oxide Calr 03 has been classified as a Mott-type insula¬ 
tor. Based on a systematic density functional theory (DFT) study with local, semilocal, and hybrid exchange- 
correlation functionals, we reveal that the Ir t 2 g states exhibit large splittings and one-dimensional electronic 
states along the c axis due to a tetragonal crystal field. Our hybrid DFT calculation adequately describes the 
antiferromagnetic (AFM) order along the c direction via a superexchange interaction between Ir‘*+ spins. Fur¬ 
thermore, the spin-orbit coupling (SOC) hybridizes the t 2 g states to open an insulating gap. These results 
indicate that Calr 03 can be represented as a spin-orbit Slater insulator, driven by the interplay between a long- 
range AFM order and the SOC. Such a Slater mechanism for the gap formation is also demonstrated by the DFT 
-I- dynamical mean field theory calculation, where the metal-insulator transition and the paramagnetic to AFM 
phase transition are concomitant with each other. 

PACS numbers: 71.20.Be, 71.70.Ej, 75.10.Lp, 71.15.Mb 


One of the most important phenomena in condensed mat¬ 
ter physics is the Mott transition driven by electron-electron 
correlations Qi. In 3d transition-metal oxides (TMOs), 
the localized 3d orbitals are responsible for the strong on-site 
Coulomb repulsion (U), leading to a Mott-Hubbard insulator 
where U splits a half-filled band into lower and upper Hub¬ 
bard bands. Surprisingly, despite weaker U in 5d TMOs due 
to the very delocalized 5d orbitals, a series of Ir oxides such as 
Sr2lr04 M, NajIrOs MM, and CalrOa |[IM1 includ¬ 
ing Ir^+ ions with five valence electrons exhibit an insulat¬ 
ing ground state. For this unusual insulating behavior of the 
5d iridates, it was proposed that spin-orbit coupling (SOC) 
splits the Ir t 2 g states into completely filled y'eff = 3/2 bands 
and a narrow half-filled y'eff = 1/2 band at the Fermi level 
(Ep), and the latter band is further split into two Hubbard sub¬ 
bands by moderate Coulomb repulsion ffl. Such a Jeff = 
1/2 Mott-Hubbard scenario has, however, been challenged by 
an alternative scenario of Slater mechanism S based on the 
single-particle band picture, where the opening of insulating 
gap in 5d TMOs is driven by a long-range magnetic order¬ 
ing |110[BIIi]. 


Here we focus on the post-perovskite CalrOs with a highly 
anisotropic geometry where IrOg octahedra share corners 
along the c axis and have common edges along the a axis (see 
Fig. 1). Recently, the nature of the ground state in Calr 03 has 
been an object of hot debate ll iMlSll . On the basis of reso¬ 
nant x-r^ magnetic scattering (RMXS) experiment, Ohgushi 
et al. il2r claimed the robustness of the y'eff = 1/2 ground 
state against structural distortions. However, a resonant in¬ 
elastic x-ray scattering (RIXS) experiment of Sala et al. lll3n 
concluded that CalrOs is not a y'eff =1/2 iridate by showing 
that the y'eff = 1/2 state is severely altered by a large tetrag¬ 
onal crystal field splitting, therefore proposing the existence 
of Mott insulator beyond the y'eff = 1/2 ground state. On the 


theoretical side, using density-functional theory (DFT) calcu¬ 
lations within the local density approximation (LDA) includ¬ 
ing SOC, Subedi 0141] interpreted the t 2 g states in terms of the 
jeS = 1/2 and y'eff = 3/2 states, and further showed that the in¬ 
troduction of U splits the y'eff = 1/2 bands into fully occupied 
lower and unoccupied upper Hubbard bands, thereby support¬ 
ing the Mott-Hubbard scenario. Contrasting with the y'eff = 
1/2 ground state generated by equally hybridizing the three 
dxy, dy^, and d^x orbitals |12|]. ab intio wave-function quantum 
chemical calculation illSll predicted the highly uneven admix¬ 
ture of xy, yz, and zx characters for the relativistic t 2 g states. 
Therefore, previous experimental and theoretical studies of 
CalrOs have not reached a consensus on the presence of the 
Jeff = 1/2 ground state, but concluded in the same way that 
CalrOs belongs to a Mott-Hubbard insulator. Despite such a 
contradiction for the nature of the ground state, it is well es¬ 
tablished ifl^UIsI] that the insulating ground state of Calr 03 
exhibits the stripe-type magnetic order with a strong antifer¬ 
romagnetic (AFM) coupling along the c axis and a weak fer¬ 
romagnetic one along the a axis (hereafter designated as the 
AFM structure). 

In this Letter, we investigate the nature of the ground state 
of Calr 03 by using comprehensive DFT calculations with lo¬ 
cal, semilocal, and hybrid exchange-correlation functionals as 
well as by including dynamical mean field theory (DMFT). 
We find that the t 2 g states are significantly split by a com¬ 
pression of IrOe octahedra along the c axis and particularly 
two t 2 g states (designated as and t^) have dominant dy. 
and characters with large band dispersions, indicating one¬ 
dimensional (ID) electronic states along the c axis. Our hy¬ 
brid DFT calculation adequately describes the delocalized t^g 
and states to stabilize the AFM order along the c axis via 


superexchange interaction between Ir^+ spins. Moreover, the 

and states with other tzg 


SOC is found to hybridize the t^g 
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FIG. 1: (Color online) Crystal structure of Calr 03 projected on the 
be plane (left side) and the ab plane (right side), a, b, and c denote 
unit vectors of the conventional unit cell. The primitive unit cell is 
indicated by the dashed lines. The large, medium, and small circles 
represent Ca, Ir, and O atoms, respectively. The bond lengths (di and 
1 ^ 2 ) and the tilt angle (6,) are indicated for Table IIS. The reference 
frame xyz is drawn in an IrOs octahedron. 


states to open an insulating gap. It is thus elucidated that the 
gap formation in CalrOs is driven by the interplay between a 
long-range AFM order and the SOC, representing a spin-orbit 
Slater insulator. This Slater mechanism for the gap formation 
in CalrOa is further supported by the DFTh-DMFT calcula¬ 
tion. 

We first study the ground state of CalrOs using the LDA, 
GGA, and hybrid DFT calculations ifl^ without SOC. Inter¬ 
estingly, the different ground states are predicted depending 
on the employed exchange-correlation functionals: i.e., the 
CA functional predicts the nonmagnetic (NM) ground state, 
while the PBE and HSE functionals favor the AEM structure 
over the NM structure by AEnm-afm = 21.1 and 85.8 meV 
per primitive unit cell, respectively. The calculated structural 
parameters of Ir-O bond lengths {d\ and c /2 in Eig. 1) and tilt 
angle (0, in Eig. 1) change little depending on the employed 
exchange-correlation functionals. Our values of c/i/c /2 ~ 0.96 
and 0, Ri 23° indicate that the IrOg octahedra are compressed 
along the c axis with a large tilt, in good agreement with an 
XRD analysis ll2^ and a previous LDA calculation 114 ] (see 
Table IIS of the Supplemental Material |21]). Eigure 2(a) 
shows the LDA band structure of the NM structure together 
with c/-orbitals projected bands and the charge characters of 
the t 2 g states. There are six t 2 g bands originating from two 
different Ir atoms within the primitive unit cell, which are 
grouped into two doublets (t^g and t^g) and two singlets (t^g 
and fjg)- We find the presence of partially occupied t 2 g states 
crossing Ep, indicating a metallic feature. It is seen in Eig. 
2(a) that (i) the tfg and doublets are significantly sepa¬ 
rated by ~ 1.5 eV at the F point and (ii) the low(high)-lying 
^2g ^^2g^ States involve dominantly dxy (dyg and dgjc) compo- 


nent(s). These different orbital characters of t^g and tP^ 


‘2g 


can 


be attributed to a tetragonal crystal field due to the compressed 
IrOg octahedra along the c axis. In addition, the t^g and t^g 
states also have dominant dyg and dgx characters with a large 


band width of ^2.5 eV, indicating ID electronic states along 
the c axis. As shown in the inset in Fig. 2(a), the charge char¬ 
acter of fj] (or t 2 g) shows a strong hybridization with 0 2/7 
orbitals along the c axis, reflecting its ID band feature. This 
ID features of the t^g and states play an important role in 
determining the electronic and magnetic properties of CalrOs, 
as discussed below. 


Figure 2(b) and 2(c) show the PBE and HSE band structures 
of the AEM structure, respectively. Here, the overall band dis¬ 
persion of the AEM structure is similar to that [Eig. 2(a)] 
of the NM structure. However, there are two conspicuous 
changes in the AEM band structure: i.e., one is a gap opening 
Eg^^^^ of the t 2 g and t^g singlets at the T point and the other 
is some separation of the two bands in the doublet. We ob¬ 
tain = 0.23 eV [Eig. 2(b)] and 0.67 eV [Eig. 2(c)] from 

PBE and HSE, respectively. To understand the microscopic 
mechanism for these changes in the AEM band structure, we 
plot the spin-polarized fi-orbitals projected bands (see Eig. 5S 
of the Supplemental Material J^ ). We find that the t^g and 


t 2 g states (or the two states in tfg) with the same spin direc¬ 
tion, localized at two different Ir sites, hybridize with each 
other, yielding an energy gain from the exchange kinetic en¬ 


ergy and the formation of Ep Il29t] . This so-called superex¬ 
change interaction | ^ 3J. ] between Ir^+ spins results in the 
stabilization of the AEM order along the c axis. Our find¬ 
ing that the t^g, t^g, and t^g states with dyg and dgx characters 
are associated with the AFM order does not support the in¬ 
terpretation of previous RMXS experiment 11211 that the y'eff 
= 1/2 state stabilizes the AFM order. In Table I, we list the 
calculated magnetic moments of Ir and O atoms in the AFM 
structure. Our HSE calculation gives a spin magnetic moment 
ms = 0.48 (0.07) for Ir (Oi) atoms, larger than mg = 0.31 
(0.05) pb computed using PBE. We note that the LDA and 
GGA tend to stabilize artificially delocalized electronic states 
due to their inherent self-interaction error because delocaliza¬ 
tion reduces the spurious self-repulsion of electrons 113211330 . 
This aspect of LDA and GGA may account for why our LDA 
calculation predicts a metallic NM ground state and our PBE 
values of AEnm-afm and mg are relatively smaller than the 
corresponding HSE ones (see Table I). 


Next, we examine the effect of SOC on the geometry and 
band structure of the AFM structure using PBE-fSOC and 
HSE-fSOC ifl^ . As shown in Table IIS of the Supplemen¬ 
tal Material, the inclusion of SOC changes little dild 2 and 0, 
by less that 0.01 and 1°, respectively. Figures 3(a) and 3(b) 
show the PBE-fSOC and HSE-fSOC band structures of the 
AFM structure, respectively. Compared to the PBE [Fig. 2(b)] 
and HSE [Fig. 2(c)] results, ^/-orbitals projected bands clearly 
show that the inclusion of SOC does not change the dxy, dyg, 
and d^ characters of the t^g, t^g and states, but gives 
rise to their strong hybridizations leading to well-separated 
t 2 g bands with hybridization gaps. This result doesn’t support 
the /eff = 1/2 ground-state picture proposed by RMXS exper¬ 
iment d, where SOC can dramatically affect the t 2 g states 
to produce the jeff = 1/2 state by equally mixing up the dxy. 
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FIG. 2: (Color online) Band structures of (a) the NM structure obtained using LDA, (b) the AFM structure obtained using PBE, and (c) 
the AFM structure obtained using HSE. The band dispersions are plotted along the path r(0,0,0) —t T(0,0.5,0) —t 7(0,0.5,0.5) —t r(0,0,0) 
—t 5(0.25,0.25,0) I?(0.25,0.25,0.5) —t r(0,0,0) —t Z(0,0,0.5) in the Brillouin zone of the unit cell [see the inset in (b)]. The energy zero 

represents the Fermi level. The bands projected onto dxy, dy^, and d„ orbitals are also displayed. Here, the radii of circle and semicircles 
are proportional to the weights of corresponding orbitals. In (a), the charge characters of and t!^g states at the T point are shown with an 
isosurface of 0.02 e/A^. 


dy^, and orbitals. We note that for the t 2 g state near Ep, 
the partial densities of states (PDOS) of the dxy, dy^, and d^x 
orbitals obtained using HSE give a ratio of 0%:50%:50%, but 
change into 8%:46%:46% for the HSE+SOC calculation [see 
Eig. 6S of the Supplemental Material]. The latter PDOS ratio 
of dxy, dy^, and d,^ is very consistent with the interpretation of 
RIXS experiment 11131] that the t 2 g state nearE/r was estimated 
as =F0.32|vy,=F > + 0.67 ((|yz,± > ^i\zx,± >). As a con¬ 
sequence of the SOC-induced hybridization, PBEh-SOC still 
does not open an insulating gap [Pig. 3(a)], but HSEh-SOC 
opens Eg = 0.32 e V [ Pig. 3(b)], close to the experimental 
data of ^0.34 eV ll 111 . Especially, the DOS [see Pig. 3(b)] 
obtained using HSEh-SOC explains well the two d-d inter¬ 
band transitions observed in RIXS experiment 11341] . It is note¬ 
worthy that the gap opening in HSEh-SOC can be realized by 
the strong interplay between a long-range APM order and the 
SOC: i.e., the enhanced separation of the t^g and bands 
in the APM structure [see Pig. 2(c)] and the SOC-induced 
hybridization are combined to yield the opening of an insulat¬ 
ing gap. Thus, we can say that the single-particle band the¬ 
ory within the HSEh-SOC scheme predicts well the magnetic- 
insulating ground state of CalrOs, being represented as a spin- 
orbit Slater insulator. 

Pigure 3(c) shows the sum of the spin and orbital mag¬ 
netic moments of Ir and O atoms, obtained using HSEh-SOC. 
The total magnetic moment for Ir, Oi, and O 2 atoms is (0.00, 
-0.14, ±0.78), (0.00, -0.01, ±0.05), and (0.00, -0.02,0.00) 
in units of /jb, respectively, showing an antiparallel alignment 
of magnetic moments along the c axis. Here, the Ir magnetic 
moments are canted along the b axis with ^11°, comparable 
with those (2-4°) reported from RMXS |[I1 and RIXS |fH 
experiments. 

Pinally, in order to confirm the Slater mechanism for the 
gap formation in CalrOs, we perform the DPTh-DMPT calcu¬ 
lation 1^ 35] for the APM and paramagnetic phases. Pigure 
4(a) shows the DPTh-DMPT results for the one-particle spec¬ 
tra. It is seen that the APM insulating phase is transformed 


TABLE I: Spin magnetic moments (in units of jAg) of Ir, Oi, 
and O 2 atoms obtained using the PBE, HSE, PBEh-SOC, and 
HSEh-SOC calculations. The calculated orbital magnetic moments 
using PBEh-SOC and HSEh-SOC are given in parentheses. Oi is 
placed in the ab plane while O 2 represents the corner-sharing O atom 
along the c axis (see Eig. 1). 



Ir 

Oi 

O 2 

PBE 

0.31 

0.05 

0.00 

HSE 

0.48 

0.07 

0.00 

PBEh-SOC 

0.16(0.08) 

0.02(0.00) 

0.00(0.00) 

HSEh-SOC 

0.45(0.34) 

0.05(0.01) 

0.01(0.01) 


into a paramagnetic metallic phase IH] around 230 K. Below 
this metal-insulator transition (MIT) temperature, the mag¬ 
netic moment of Ir atom sharply increases [see Pig. 4(b)], 
indicating that the MIT and the paramagnetic to APM phase 
transition are concomitant with each other. The calculated 
spectral function for the APM ground state is displayed in Pig. 
7S of the Supplemental Material im, together with the super¬ 
imposition of the HSEh-SOC band structure. Both results are 
in good agreement with each other. Therefore, both the com¬ 
prehensive DPT calculations and the DPTh-DMPT calculation 
reach a consistent conclusion that the gap opening is induced 
by the APM order, representing a Slater-type insulator. 

In conclusion, our comprehensive DPT calculations with 
local, semilocal, and hybrid exchange-correlation functionals 
clarihed the effects of tetragonal crystal held, APM order, and 
SOC on the t 2 g states of CalrOs. We found the large tetragonal 
crystal held splitting of the t 2 g states, resulting in the forma¬ 
tions of t 2 g doublet (with and d^x characters), t^g doublet 
(with dxy character), and t^g and singlets (with dy^ and d^x 
characters). We also found that the insulating-gap opening 
is formed by the interplay between a long-range APM order 
and the SOC, representing a spin-orbit Slater insulator. This 
single-particle based Slater picture is further supported by the 
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FIG. 3: (Color online) Band structure of the AFM structure, obtained using the (a) PBE+SOC and (b) HSE+SOC calculations. The DOS 
obtained using HSE+SOC is also given in (b). The theoretical estimation of the B and C features observed from RIXS experiment Q 
is drawn in DOS. In (c), the sum of the spin and orbital magnetic moments obtained using the HSE+SOC calculation is drawn with three 
components (my, mi,, me) along the a, b, and c axes. Here, mi is calculated by integrating the corresponding component of magnetic moment 
inside the PAW sphere with a radius of 1.4 (0.8) A for Ir (O). 
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FIG. 4: (Color online) (a) Density of states with respect to temper¬ 
ature, obtained using the DFT+DMFT calculation. The calculated 
magnetic moment (/r) of Ir atom as a function of temperature is plot¬ 
ted in (b), together with spin (/y^) and orbital (^ii) components. 


DFT+DMFT calculation. Our findings are anticipated to stim¬ 
ulate further experimental and theoretical studies for other iri- 
dates such as Sr 2 lr 04 and Na 2 lr 03 lIsUI^ which have 
been proposed to be a y'eff = 1/2 Mott-Hubbard insulator. 
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1. Calculated structural parameters in the optimized structure of CalrOs 


TABLE IIS: Ir-O bond lengths (di and d 2 in Fig. I) and tilt angle (0, in Fig. 1) obtained using the LDA, PBE, and PBE+SOC calculations, in 
comparison with the experimental data. For the HSE and HSE+SOC calculations, we used the optimized structures using PBE and PBE+SOC, 
respectively. 



di (A) 

diik) 

e,(°) 

EDA 

1.969 

2.035 

22.2 

PBE 

1.973 

2.049 

22.5 

PBE+SOC 

1.976 

2.052 

22.6 

Experiment (Ref. [20]) 

1.972 

2.049 

22.3 


2. Spin-polarized d-orbitals projected band structures 

To understand the microscopic mechanism for the antiferromagnetic (AFM) spin ordering in CalrOs, we plot the spin- 
polarized tf-orbitals projected band structures obtained using the HSE calculation in Fig. 5S. It is seen that the doubly degenerate 
(spin-up and spin-down) states are localized at Iri and Ir 2 sites. Since electronic states with the same spin direction can hybridize 
with each other, the hybridization takes place between the spin-up (or spin-down) states localized at the Iri and Ir 2 sites, yielding 
not only an energy gain from the exchange kinetic energy but also the formation of hybridization gap (E|'^'^^) or the enhanced 
energy splitting. This so-called superexchange interaction between Ir^+ spins results in the stabilization of the AFM order along 
the c axis. 
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FIG. 5S: Spin-polarized d-orbitals projected band structures of the AFM structure obtained using the HSE calculation. The bands projected 
onto dxy, dyz, and d^ orbitals are displayed with circle and semicircles whose the radii are proportional to the weights of corresponding orbitals. 
Iri Ir 2 are indicated in the Fig. 1 of the main text. The energy zero represents the Fermi level. 




























































7 


3. Partial densities of states of the dxy, dy^, and d^ orbitals 

We plot the partial densities of states (PDOS) of the dxy, dyz, and d^x orbitals using the HSE and HSE+SOC calculation in Eig. 
6S. Eor the t 2 g state above Ep, the PDOS obtained using HSE give a ratio of 0%:50%:50%, but they change into 8%:46%:46% 
for the HSE+SOC calculation. The latter PDOS ratio of the dxy, dy^, and d^x orbitals is very consistent with the interpretation of 
RIXS experiment [13] that the t 2 g state near Ep was estimated as =F0.32|.ry, =F > + 0.67 ((|yz,± > ^i\zx,±. >). 
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FIG. 6S: PDOS of the dxy, dy^, and dzx orbitals obtained by the (a) HSE and (b) HSE+SOC calculation. 


4. Computational details of the DFT+DMFT calculation 

The dynamical mean field theory (DMET) calculations are performed in a fully charge self-consistent way by combining 
with density functional theory (DPT) using the projection-embedding implementation [22, 35]. The DPT calculations within 
the generalized gradient approximation (GGA) [23] are carried out using the full-potential linearized augmented plane-wave 
(PLAPW) method as implemented in the WIEN2K package [24]. In the DMET part, we solve the quantum impurity model 
using the continuous-time quantum Monte Carlo method [25, 26]. The local correlated orbitals are projected from a 20eV 
energy window around the Permi level. The double-counting functional is expressed as ‘Pdc{nd) = ^dEdc, where Edc = U{rij — 
1 /2) — y/2(«® — 1) and the nominal occupancy n® of Ir ion in CalrOs is 5. Due to the large energy window, the local Coulomb 
interaction parameters are relatively system-independent and estimated [27] asU ^ 4.5 eV and J Ri 0.8 eV [28]. We set different 
local coordinates on each Ir atom in order to simulate the antiferromagnetic structure. With the crystal-field and spin-orbit 
coupling interactions considered, the effective 7eff =1/2 states are constructed as, 

\dxy,To) + '^{\dyg,±a)Ti\dxz,±o)) (1) 

where a labels the spin [28]. The value of y at low-frequency limit is adopted in our transformation in which y(® = 0) r; 1.16. 
Then the magnitudes of magnetic moments along the magnetization axis can be calculated as, {/jj) = I'f'An /3 and (pj) = 
( 4 y 2 _ 3)A«/3, where An is the difference of the occupation numbers of |'|t±i/ 2 ) [28]. 
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5. Spectral function of the AFM state obtained using DFT+DMFT 

The calculated spectral function of the AFM ground state is displayed in Fig. 7S, together with the superimposition of the 
HSE+SOC band structure. Both results are in good agreement with each other. 



FIG. 7S: Spectral function of the AFM state obtained using DFT+DMFT. Dotted lines denote the band structure obtained using FISE+SOC. 
Here, the DFT+DMFT spectral function and the HSE+SOC band structure are aligned to have the same position at the topmost valence state 
atF. 









